Low-temperature triple-alpha rate in a full three-body model 
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A new three-body method is used to compute the rate of the triple-alpha capture reaction which is 
the primary source of 12 C in stars. In this work, we combine the Faddeev hyperspherical harmonics 
and the R-matrix method to obtain a full solution to the three-body a + a + a continuum. Particular 
attention is paid to the long range effects caused by the pairwise Coulomb interactions. The new rate 
agrees with the NACRE rate for temperatures greater than 0.07 GK, but a large enhancement at 
lower temperature is found (~ 10 12 at 0.02 GK). Our results are compared to previous calculations 
where additional approximations were made. We show that the new rate does not significantly 
change the evolution of stars around one solar mass. In particular, such stars still undergo a red- 
giant phase consistent with observations, and no significant differences are found in the final white 
dwarfs. 
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The triple-alpha reaction is key to the production of 
elements with mass A > 12. It plays a role in many as- 
trophysical processes, from stellar evolution to explosive 
scenarios. Since the early days when Fred Hoyle postu- 
lated that there should be a + resonance close to the 
3a threshold to justify the observed abundances of 12 C 
in stars pQ, a state experimentally confirmed by Dunbar 
et al. [2], the accepted mechanism for the triple-alpha 
capture has been sequential: a + a —> 8 Be(0^) followed 
by 8 Be(0^) + a ->■ 12 C(0.[). However, if the energy of the 
alpha particles in the stellar environment is insufficient 
to populate the narrow resonances, a direct three-body 
capture becomes the favored path. The most recent com- 
pilation of charge-particle induced thermonuclear rates 
(NACRE) [3 extrapolates the sequential model to low 
energy, an approach that should be questioned j3] [5] . Re- 
cent theoretical works [JJ |S] predict enhancements in the 
triple-alpha rate by 7 to 20 orders of magnitude when 
compared to NACRE, for temperatures around 0.02 GK. 
Because the triple-alpha rate, even at these low temper- 
atures, can have dramatic and widespread effects in as- 
trophysics, it is critical to resolve the large discrepancies 
found between [3], [3] and [5]. This is the aim of the 
present work. 

The Hoyle resonance dominates the triple-alpha rate 
at T ~ 0.1 GK, while other resonances play a role at 
T ~ 1 GK relevant in explosive environments. At these 
high temperatures, a 12 C(0^) state around 3 MeV above 
the 3a threshold [5] contributes significantly to the rate 
but the main path for the capture is through a high-lying 
12 C(2^) resonance, which decays directly to the ground 
state of 12 C [7|. Uncertainties in the properties of this 2 ~\ 
state have been recently addressed [8] . It is important to 
note that this 2 + resonance has little influence at very 
low temperatures. Even if the triple- alpha resonant rate 



remains more uncertain than desired [5], it is far better 
understood that the non-resonant counterpart. 

At temperatures below T w 0.07 GK, the 3 a-particles, 
having no access to intermediate resonances, fuse di- 
rectly. The description of this process requires a good 
solution to the three-body scattering equation including 
Coulomb effects. Recent work by the Kyushu group [3] 
makes use of the continuum discretized coupled channel 
(CDCC) method [TU] which expands the full three-body 
wavefunction in terms of the continuum states in the two- 
body subsystem (in this case 8 Be). In this method, reso- 
nant and non-resonant processes are treated on the same 
footing. The CDCC method has been widely used in 
nuclear reactions but the application to the triple-alpha 
capture is particularly challenging because the charged 
particle reaction occurs at large distances and is domi- 
nated by the Coulomb interactions. The CDCC results 
[4], when compared to NACRE, show a 20 order of mag- 
nitude increase in the reaction rate at T = 0.02 GK. Soon 
after these results became available, numerous astrophys- 
ical studies were performed to understand the implica- 
tions of such a dramatic change. Stellar evolution cal- 
culations [HI [T2] demonstrate that the CDCC rate sup- 
presses He flashes in the asymptotic giant branch, and 
can result in the disappearance of the red giant phase 
in low and intermediate mass stars. Studies of He igni- 
tion in accreting white dwarfs |13j and accreting neutron 
stars O [T5] show that the CDCC rate is barely consis- 
tent with observations of Type la supernovae and type I 
X-ray bursts, respectively. 

The CDCC calculations in [3] represent the first at- 
tempt to take into account the non-resonant contribu- 
tion to the triple alpha cross section at very low relative 
energies, where it is expected to dominate. The main 
drawback of the CDCC method is that, in general, the 



truncated CDCC wavefunction does not have the correct 
asymptotic form for the scattering of three charged par- 
ticles [IB] . It is known that the asymptotic form is nat- 
urally expressed in terms of hyperspherical coordinates 
[IB] , Recent attempts to solve the a + a + a problem 
by the Madrid- Aarhus collaboration [T7] [T5] employ the 
hyperspherical adiabatic expansion method. Limitations 
in the numerical techniques did not allow for determin- 
ing the rate below T 0.1 GK. In this work we discuss 
another method based on the hyperspherical coordinates 
which does provide improved accuracy and enables us to 
reach temperatures down to 0.01 GK. 

An alternative, rather schematic, approach was re- 
cently proposed [5]. Inspired by the low-energy extrapo- 
lation of the sequential process, implemented in NACRE, 
Ref. [5j proposed an extrapolation of the three-body 
Breit-Wigner form for the contribution to the three-body 
capture that does not go through the 8 Be resonance. 
This method is here referred as BW(3B). Results for the 
triple-alpha reaction rate in [5] show a modification of 
the temperature dependence of the rate at low tempera- 
ture amounting to 7 orders of magnitude increase of the 
triple-alpha rate at 0.02 GK, when compared to NACRE. 

The works mentioned before [U [18] approach the 
triple-alpha reaction as a three-body problem. Presently, 
a fully microscopic approach to the non-resonant triple- 
alpha capture is not feasible. In this work we also follow 
the three-body approach, and expect microscopic cor- 
rections to this approach to be small compared to the 
three-body effects discussed in [JJ [5] . 

The reaction rate R a b c for a + b + c — > D + 7 at a 
given energy E can be related to the cross section for 
photo-dissociation <r 7 through: 
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with E 1 the photon energy, E the three-body energy, /n a b 
and Hab,c the reduced masses of the a + b and (ab) + c 
systems, respectively, p the number of identical particles, 
Na Avogadro's number and G a b c ,D a statistical factor 
depending on the spins of the nuclei a,b,c and D [S]. 
The energies are related by E 1 = E + \Er>\, where Ed 
is the energy of the bound final state of D relative to 
the 3a threshold. The relevant astrophysical quantity is 
the energy averaged reaction rate (R a bc)(T), obtained by 
integrating R a b c (E) over the Maxwell-Boltzmann distri- 
bution. 

At low energies, the photo-dissociation process is usu- 
ally dominated by a single transition of multipolarity A, 
so 
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depending on the corresponding electromagnetic tran- 
sition strength dB(EX)/dE. For the triple-alpha reac- 
tion at low energies, the reaction proceeds through a 
quadrupole transition from the + continuum to the 2f 
bound state in 12 C. 

In order to calculate the transition strength 
dB(E2)/dE at low relative energies, we need to 
solve the bound three-body problem to obtain the 
2f state and to solve the full scattering three-body 
problem for the + triple-alpha continuum. To solve the 
three-body Faddeev equation, we use the hyperspherical 
harmonic (HH) method |22j . Hyperspherical coordi- 
nates (hyper-radius p = \J xf + and hyper-angle 
9; = tan—) are introduced as functions of the scaled 
Jacobi coordinates Xi,yi, with i = 1,2,3 22J. The 
HH method expands the three-body wavefunction as 
sums of products of hyper-radial functions x-y{p) an d 
angular functions fl 7 (9,x,y) where 7 = {£, Ki, l Xi , l yi }. 
This expansion introduces a new quantum number, the 
hyper-momentum Ki, and sufficiently many Ki values 
need to be used to obtain convergence. The orbital 
angular momenta l Xi , l Vi refer to the Jacobi coordinates 
Xi,yi. Introducing the HH expansion in the Faddeev 
equations gives the set of coupled channels equations in 
the hyper-radius coordinate: 
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where Aj = Ki + 3/2 and m is a scaling mass. The 
coupling potentials are defined as the sum of three pair- 
wise interactions Vjk{p, 0) plus a three-body force V^bip), 
integrated over all variables but p: 
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We take the same interactions for the three-body Hamil- 
tonian as in [17] . After minor adjustments, this repro- 
duces the 8 Be ground state, and the 12 C 2+ bound state 
and O2 Hoyle resonance. We explicitly take into account 
the symmetrization of the system, thereby reducing the 
set of coupled equations. 

The solution for the 2f bound state is obtained by solv- 
ing Eq. ([3| for negative energy E. The numerical method 
is implemented in FACE [23], which is specifically de- 
signed for three-body bound states. The + continuum 
states are found by using the R-matrix method in the 
HH basis (HHR), originally developed in the context of 
core+n+n scattering [24j . a problem with different sym- 
metries and no Coulomb interactions. The HHR method 
generates an orthonormal basis set to solve Eq. ([3| in 
a box with hyper-radius p m , by fixing the logarithmic 
derivative there. We expand \i m terms of the HHR 
basis. Due to the long range Coulomb forces, numeri- 
cal solutions are needed out to the very large distances 
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(p 100 fm) where it becomes safe to ignore off-diagonal 
couplings and matching to known Coulomb functions can 
be performed [16]. Extending the orthonormal basis to 
the required large p would be impractical, so we use an R- 
matrix propagation technique |25j out to an asymptotic 
Pa Pm- We divide the interval from p m to p a into 
sectors labelled p, the sizes of which are sufficiently small 
that the interactions within can be taken as constant. 
The propagating functions G p are then given in terms of 
complex exponentials |25j . and used to propagate from 
sector p—1 to the next: 

[_r p_i - c?{piy L )]-^ . (5) 

This gives the R-matrix on the right side of the bound- 
ary BP from B p ~ l on the left. The coordinate sub- 
scripts p/f, pl imply evaluations on the right and the 
left side of the boundary respectively. While without 
Coulomb this propagation method was numerically sta- 
ble, including the Coulomb interaction introduced nu- 
merical instabilities in the three-body scattering wave- 
function, in the turning point region. To overcome these 
instabilities, we introduced screening in the off-diagonal 
potentials Eq.Q using a Woods-Saxon multiplying fac- 
tor [l + exp(( j o-p scrcon )/a scroon )] _1 . Our calculations use 
a sufficiently large screening radius such that dB(E2) / dE 
for E = 0.05-0.5 MeV becomes independent of p scre en 

and ^screen- 

Convergence of dB(E2) /dE was checked for all the 
relevant parameters. At relative energies ps 0.1 MeV, 
dB{E2) / dE picks up contributions from the wave func- 
tions at hyper-radii p = 5-30 fm. However, in order 
to get a reliable three-body scattering wavefunction, the 
matching to the asymptotic form needs to be performed 
at p > 

pscrccn- 

The results here presented include up 
to K — 26 in the HH expansion, 50 poles in the R- 
matrix expansion, an R-matrix box size of 50 fm and 
R-matrix propagation out to 3000 fm. In addition we use 
Pscrocn = 800 fm and a scree n = 10 fm to screen the off- 
diagonal couplings. The only limiting factor for full con- 
vergence is the truncation in K max . We carefully studied 
the exponential behavior of dB(E2) / 'dE as a function of 
ifmax, for the relative energy range E = 0.01-0.15 MeV 
(the relevant region for rates below 0.07 GK). We con- 
clude that our results with K max = 26 differ from those 
obtained from the extrapolation of K max — > oo by a fac- 
tor of at most 2. Our energy averaged rates are plotted 
in FigJT] as a function of temperature (solid line). The 
factor of 2 error on the rate at low energy cannot be seen 
given the scales in Fig[T] 

Our method to determine the three-body scattering 
state can be checked for the pure Coulomb case. In that 
case, no resonances are present in either the two-body or 
three-body systems and the capture is necessarily direct. 
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FIG. 1: (Color online) Different evaluations of the triple- 
alpha reaction rate: comparing the Hyperspherical Harmonic 
R-matrix method (solid) with NACRE (dotted), CDCC 
(dashed) and the three-body Breit Wigner (dot-dashed). 



For the pure Coulomb case, the propagation technique 
is stable and we were able to unambiguously determine 
that no error is introduced by screening the off-diagonal 
couplings. In addition, there is an analytic solution of 
Eq. ^ for the pure diagonal Coulomb couplings which 
we use to test our implementation. The rate for diagonal 
Coulomb differs from the full rate by 10 orders of mag- 
nitude, thus demonstrating the importance of including 
the Coulomb effects correctly. The off-diagonal Coulomb 
couplings are important, but their major contribution 
to the rates in the range T = 0.01-1.00 GK occurs for 
p < 100 fm. 

In Fig. [I] we compare our results (solid) with those 
obtained previously. The rate obtained from the full so- 
lution of the HHR equations agree with NACRE (dotted 
line) above 0.07 GK. Although there is a slight reduction 
of the rate below 0.07 GK, we find a pronounced increase 
of the rate for T < 0.06 GK, accompanied by a very 
different temperature dependence. In some ways, the re- 
sults obtained assuming a rate extrapolated with a three- 
body Breit- Wigner form BW(3B) (dot-dashed line) [5] 
are qualitatively similar, but the BW(3B) treatment en- 
hances the rate to a smaller degree. Although we do find 
a strong enhancement at low energy, it is not as strong as 
that seen in the CDCC result of [4] (long-dashed). While 
the deviation of the CDCC results from NACRE is seen 
all the way up to 0.11 GK, the HHR result agrees with 
NACRE for T > 0.07 GK. This aspect turns out to be 
crucial in the repercussions for astrophysics. 

The kink observed (Fig. [I]) in the full rate around 
T ps 0.06 GK is a signature of two processes with dif- 
ferent temperature dependence. Above T pa 0.06 GK 
the process through the Holye state dominates, while 
below there is mostly non-sequential (direct) capture. 
The calculations in [5] exhibit this same feature. The 
temperature dependence is best illustrated by evaluating 
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d\n(R aaa ) /dbxT, as shown in Table[T] For T < 0.06 GK, 
the HHR rate has a much weaker temperature depen- 
dence than the NACRE rate, as is evident in Fig. [T] For 
T > 0.06 GK, the HHR rate follows that of NACRE, and 
their temperature sensitivities are basically equal (Ta- 
ble§. 

We have tested the sensitivity of the rate to the nuclear 
interaction and the three-body force by using the same 
Hamiltonian as [J] to calculate the + continuum. The 
introduction of a three-body interaction was necessary to 
reproduce the relevant Hoyle state. This modified Hamil- 
tonian increases the rates at T = 0.02 GK by 4 order 
of magnitude, compared to those obtained with the Ali- 
Bodmer interaction }26| . However, while the Ali-Bodmer 
interaction (used in [17] and also here), reproduces the 
a-a phase shifts, that of [4] does not. 

Next we study the effect of the strong enhancement of 
the triple-alpha rate at low temperature on stellar evo- 
lution. We use the mesa (Modules for Experiments in 
Stellar Astrophysics) code library [37], but replace the 
NACRE rate by the HHR rate for T < 0.1 GK. We 
explore the evolution of progenitors with mass M = 
0.8 M Q , 1.0 M Q , and 1.25 Mq, all with a standard so- 
lar composition, using the setups described in Ref. [57J 
§ 7.1]. In general, He burning occurs in low-mass stars 
only for T > 0.07 GK, and therefore the changes to the 
stars are minimal. Figure [2] illustrates the post main 
sequence evolutionary track, in the luminosity-effective 
temperature plane, for a star of mass 1M Q . Tracks for 
the standard NACRE rate (dashed line) and the HHR 
rate (solid line) are shown. Both stars go through a red- 
giant phase until He ignites in the semi-degenerate core, 
unlike the case calculated with the CDCC rate, for which 
the red-giant phase vanishes [IT]. Only after cessation 
of core He burning and the onset of H and He thermal 
shell flashes do small differences appear (shown in Fig{2]): 
the enhanced HHR rate at T < 0.06 GK tends to pro- 
duce a slightly larger convective zone during helium shell 
flashes. Both rates lead to final white dwarfs with iden- 
tical masses and carbon-oxygen ratios. The evolutionary 
tracks for a 1.25 Mq star behave similarly to the 1M 
case, while the tracks for the 0.8 Mq star are identical for 
both the NACRE and HHR rates. Further astrophysical 
studies to understand the implications of the new rate in 



TABLE I: Temperature sensitivity of triple-alpha rate 
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FIG. 2: (Color online) Evolutionary track (luminosity vs. sur- 
face effective temperature) of a one solar mass star with solar 
composition, for the HHR rate (solid line) and the NACRE 
rate (dashed line). The evolution is identical for both from 
when H fuses to He in the core ("main sequence") thorough 
the formation of a degenerate He core ("giant branch") and 
the ignition of He in the core ("core He flash"). Small dif- 
ferences are seen when thermally unstable H and He burn- 
ing occurs in a shell about a degenerate C/O core ("thermal 
pulses"), but there is no difference in the final white dwarf's 
mass or composition. 

explosive scenarios will be carried out in the near future. 

In summary, we have developed a framework combin- 
ing the Hyperspherical method, the R-matrix expansion, 
and the R-matrix propagation and screening technique 
(named here HHR) to calculate the triple-alpha reaction 
rate at low temperatures, T < 0.1 GK. Our results show 
a strong enhancement of the rate below 0.06 GK, while 
reproducing the NACRE rate at higher temperatures. A 
longer publication with a detailed account of the method 
and further analysis of the reaction dynamics is under- 
way. 
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